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Abstract 
The dynamics of the plasma in the inner regions of an accretion disk around accreting millisecond 
X-ray pulsars is controlled by the magnetic field of the neutron star. The interaction between an 
accretion disk and a strong magnetic field is not well-understood, particularly at low accretion rates 
(the so-called "propeller regime" ) . This is due in part to the lack of clear observational diagnostics to 
constrain the physics of the disk-field interaction. Here we associate the strong ~ 1 Hz modulation seen 
in the accreting millisecond X-ray pulsar NGC 6440 X-2 with an instability that arises when the inner 
edge of the accretion disk is close to the corotation radius (where the stellar rotation rate matches the 
Keplerian speed in the disk). A similar modulation has previously been observed in another accreting 
millisecond X-ray pulsar (SAX J1808. 4-3658) and we suggest that the two phenomena are related and 
that this may be a common phenomenon among other magnetized systems. Detailed comparisons with 
theoretical models suggest that when the instability is observed, the interaction region between the 
disk and the field is very narrow - of the order of 1 km. Modeling further suggests that there is 
a transition region (~ 1 — 10 km) around the corotation radius where the disk-field torque changes 
sign from spin up to spin down. This is the first time that a direct observational constraint has been 
placed on the width of the disk-magnetospherc interaction region, in the frame of the trapped-disk 
instability model. 
Subject headings: X-rays: binaries — pulsars: individual (NGC 6440 X-2) 



1. INTRODUCTION 

The serendipitous discovery of the accreting millisec- 
ond X-ray pulsar (AMXP) NGC 6440 X-2 has pre- 
sented new challenges to our understanding of magnc- 
tospheric accretion. This source was discovered in 2009, 
and is, on average, the faintest known AMXP (Altami- 
rano et al. 2010). The neutron star spins at 206 Hz 
and orbits around its white dwarf companion in about 
1 hr(Altamirano et al. 2010). Between 2009 and 2011 
it underwent 10 outbursts, all of which were unusually 
short (2-3 days) with a very short recurrence time (~ 1 
month; Hcinkc et al. 2010) and peak luminosities of 
about 10 36 ergs _1 . Subsequent time-series analysis re- 
vealed a strong modulation at 1 Hz during an outburst 
(Patruno et al. 2010). This modulation is similar to that 
observed in SAX J1808. 4-3658 in the decay tail of a 
number of outbursts (van der Klis et al. 2000; Patruno 
et al. 2009a). 

Both sources show the modulation at moderately low 
luminosities, around what is known as the "propeller 
regime" . Accretion disks surrounding strong stellar mag- 
netic fields are disrupted close to the star, where the field 
is strong enough to control the gas dynamics. The loca- 
tion where the disk ends is usually termed the magne- 
tosphcric radius (r m ). If the star is spinning rapidly or 
the accretion rate is sufficiently low, the magnctospheric 
radius will lie outside the corotation radius, the point at 
which the Keplerian disk rotates at the same rate as the 
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The rapidly-spinning magnetic field then acts as a cen- 
trifugal barrier that inhibits accretion onto the star. It is 
usually assumed that the gas in the disk will then be ex- 
pelled from the system, much like a propeller (Illarionov 
& Sunyaev 1975). 

Recent new work, however, has suggested this picture 
is incomplete. Spruit & Taam (1993) and Rappaport 
et al. (2004) noted that in order for a large amount of 
mass to be expelled from the system the inner edge of 
the disk must be (at minimum) substantially outside the 
corotation radius. If not, the interaction between the 
magnetic field and disk will prevent gas from accret- 
ing but not inject enough energy to expel it from the 
disk. Furthermore, the coupling between the magnetic 
field and disk will transfer angular momentum from the 
star to the disk, which will both spin down the star and 
change the radial gas density distribution in the inner 
disk region. This disk solution, known as a "dead disk" , 
was first described by Siuniacv & Shakura 1977, but has 
been little-examined since then. 

When a dead disk forms in the inner disk regions, 
D'Angelo & Spruit (2010, 2012) demonstrated that con- 
tinued accretion from large radii will often drive an ac- 
cretion instability. The instability sets in when the inner 
edge of the disk moves outside r co , halting accretion onto 
the star. As gas continues to accumulate from the outer 
regions of the disk, the inner edge of the disk gradually 
overcomes the centrifugal barrier created by the magnetic 
field and eventually crosses r co . Once the disk inner edge 
has moved inside r co , the accumulated reservoir of gas is 
depleted and the inner edge again moves outside r co , al- 
lowing a new cycle to begin. 
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TABLE 1 
RXTE and Swift observations OF NGC 6440 X-2 



Outburst Start 


1 Hz mod. 


Instrument 


(Obs-Id) 


04-Jun-2009 


NA 


Swift/XRT 


31421001 


30-Jul-2009 


NA 


RXTE/PCA 


Bulge-Scan 


30-Aug-2009 


Y 


RXTE/PCA 


94044-04-02-00 


Ol-Oct-2009 


NA 


RXTE/PCA 


Bulge-Scan 


28-Oct-2009 


NA 


RXTE/PCA 


Bulge-Scan 


21-Mar-2010 


Y 


RXTE/PCA 


94315-01-12-00; 
94315-01-12-02 


12-Jun-2010 


Y 


RXTE/PCA 


94315-01-14-00 


04-Oct-2010 


Y 


RXTE/PCA 


94315-01-25-00 


23-Jan-2011 


Y 


RXTE/PCA 


96326-01-02-00 


06-Nov-2011 


Y 


RXTE/PCA 


96326-01-35-00 



Note. — The outbursts discussed in this paper arc highlighted 
in bold text. 

Building on earlier work by Spruit & Taam (1993) 
(which did not incorporate the dead disk solution), 
D'Angelo & Spruit (2010, 2012) studied this instabil- 
ity extensively, demonstrating that it in fact exists in 
two different accretion regimes, one at very low accre- 
tion rates, and one at higher fluxes compatible to those 
observed in some AMXPs at the onset of the propeller 
stage. 

The dead disk model also makes a radically different 
prediction for disks in quiescence from the standard one. 
In the conventional model, the inner edge of the disk 
is a straightforward function of accretion rate M, mov- 
ing outwards as M decreases. In contrast, the dead disk 
model predicts that even at very low accretion rates the 
inner edge of the disk will not move far outside r co , and 
a significant amount of gas will remain in the inner disk 
region even when M ~ 0. D'Angelo & Spruit (2011) 
demonstrated that in this case a 'trapped disk' can de- 
velop, in which the disk inner edge stays trapped close 
to r co and acts as an efficient sink for angular momen- 
tum, spinning down the accretor considerably. This disk 
solution will also likely alter the duration of individual 
outbursts, which in the ionization disk instability model 
( see e.g. Lasota 2001 for a review), is affected by both 
the surface density in the inner disk regions and the feed- 
back of radiation from accretion onto the central object 
back onto the disk. 

In this paper we present an analysis of the unusual 
outburst behaviour seen in NGC 6440 X-2 and the 1 Hz 
modulation and argue that both are manifestations of an 
unstable dead disk, following the model of D'Angelo & 
Spruit (2010). Using data from the Rossi X-ray Timing 
Explorer (RXTE) we identify the 1 Hz modulations in six 
different outbursts and we study its spectral and timing 
properties, demonstrating that these are consistent with 
coming from accretion rate modulation onto the neutron 
star. We propose that the 1 Hz modulations seen in 
NGC 6440 X-2 and SAX J1808.4-3658 have the same 
origin and are a new and unique diagnostic to study the 
disk-magnctosphcrc interaction. 

2. X-RAY OBSERVATIONS 

We used all pointed RXTE observations recorded with 
the Proportional Counter Array (PCA; Jahoda et al. 
2006) that cover a time range from August 30, 2009 
(MJD 55073) to November 18, 2011 (MJD 55883). We 



used the Standard 2 mode data (16-s time resolution) 
to calculate the X-ray lightcurve. The background 
counts are subtracted and dead-time corrections arc 
made. The energy-channel conversion is done by using 
the pca_e2c_e05v04 table provided by the RXTE team. 
In order to correct for the gain changes and the differ- 
ences in effective area between the different PCUs, the 
flux is normalized to the Crab Nebula values (see Kuulk- 
ers et al. 1994; van Straaten et al. 2003) that are closest 
in time but in the same RXTE gain epoch (Jahoda et al. 
2006). 

For the timing analysis we used all GoodXcnon and 
Event data, with a time resolution of 2 -20 and 2 -13 s, 
respectively. The GoodXenon data were rebinned to the 
same resolution of the Event data. No background sub- 
traction or dead-time correction was applied to the data 
before calculating the power spectra. The Poissonian 
noise level was determined by taking the average power 
between 3000 and 4000 Hz, a region in the power spectra 
dominated by counting statistics noise. This mean value 
was then subtracted from the power spectra. The back- 
ground was calculated for each observation (Obs-Id) with 
the FTOOL pcabackest over the entire energy range. 

We used 256s-long segments to calculate power spec- 
tra with no energy selection (i.e., using all 256 energy 
channels) . The frequency boundaries for the power spec- 
tra are therefore 1/256 Hz and 4096 Hz. The power 
spectra were normalized in the rms normalization (van 
der Klis 1995) which gives the power density in units 
of (rms/mean) 2 Hz -1 . Consistent with Patruno et al. 
(2009a), we define the fractional rms amplitude between 
v\ and V2 as: 
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P(v)i 
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and calculate the errors from the dispersion of the data 
points in the power spectra. 

To confirm the presence of coherent X-ray pulsations 
at the spin frequency (sa 206 Hz) of the neutron star, we 
barycentered the photons arrival times with the FTOOL 
f axbary using the most precise Chandra astrometric po- 
sition available (Heinke et al. 2010) and the JPL DE-405 
solar system ephemeris. We then generated dynamical 
power spectra with 256 s long data segments and over- 
lapping intervals of 64 s. As is common practice in co- 
herent timing studies of AMXPs, we here retain only the 
energy channels between 5 and 37 (absolute channels). 

Each power spectrum is inspected by eye to look for 
the presence of the 1 Hz modulation which was already 
observed and reported on June 12th 2010 (Patruno et al. 
2010). The 1 Hz modulation was observed in the power 
spectrum of SAX J1808. 4-3658 as a broad feature that 
can be described by one or multiple Gaussians (Patruno 
et al. 2009a). Therefore we follow the same procedure 
here and fit the 1 Hz modulation in the power spectrum 
with one or more Gaussian functions: 



P(y) = Acxp[(v - v Q ) 2 / C 2 



(3) 



where A is the amplitude of the modulation, vq the cen- 
troid frequency and the constant C is related to the 
FWHM via the expression: 



FWHM = 2C ■ [ln(2) 
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(4) 
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Fig. 1.— The 2-16 keV X-ray lightcurve of NGC 6440 X-2 from 
Aug 30, 2009 (MJD 55073) until November 18, 2011. The plot 
shows the occurrence of six outbursts highlighted with a number. 
The 1 Hz modulation is detected in all six outbursts with the flux 
that is always above ~ 2mCrab. 

The choice of a Gaussian distribution rather than a 
Lorentzian distribution (which is more typical in QPOs 
fitting) is made because the power spectrum shows a 
rapid decline at lower frequencies, leading to a systematic 
overestimate when fitting with a Lorentzian function. 



o 




Fig. 2. — Ten-second sa 2 — 16 keV lightcurve with the 1 Hz mod- 
ulation clearly observed. The lightcurve is created by using Event 
data (Obsld: 94315-01-14-00) to have a sufficient time resolution 
to display the 1 Hz modulation. 

We study the energy dependence of the 1 Hz mod- 
ulation by selecting the observation with the highest 
signal-to-noise (Obsld 94315-01-14-00) and splitting the 
data into four energy bands (2.47—5.31, 5.71—9.81, 
10.22-15.60 and 16.01-28.67 keV). The intrinsic faint- 
ness of NGC 6440 X-2 does not allow a splitting of the 
data into more sub-bands. We do not select any energy 
above « 28 keV as the count rate is consistent with the 
background. 

We also investigate the coherence of the time vari- 
ability in different energy bands, closely following the 
method outlined in Vaughan & Nowak (1997) and 
Scaringi et al. (2013). The coherence between two time 



series is defined as the degree of linear correlation be- 
tween two Fourier frequencies in two different lightcurves 
observed simultaneously but in different energy bands. 
Once each lightcurve is corrected for the Poissonian noise 
component, then it is possible to calculate the intrinsic 
coherence function 7j 2 „ t . In simple terms, the intrinsic 
coherence function describes how strongly the variability 
observed in one energy band is causally connected to the 
variability seen in the other band. If the the degree of 
coherence is high, then the variability seen in one energy 
band is correlated to the variability observed in the other 
band. 

In our measurements, we use again the observation 
94315-01-14-00 (w 9 ks exposure) and define two energy 
bands: 2.47-7.76 keV (soft band "x") and 8.17-28.67 
keV (hard band "y"). The two bands are defined so that 
the counting statistics is similar in both time series. We 
split the two time series in N = 33 data segments of 
length 256 s each. Each z-th time series Xi(t) and yi(t) 
(with i = 1, ...,N) is Fourier transformed into Xiiy) and 
Yi{v). We then calculate the cross-spectra between each 
pair as: 

C i (u)=X*(v)Y i (p) (5) 

and the power spectral density of each time series: 
P x (v) = \X(v)fandP y (v) = \Y(v)\ 2 

In an ideal linear system with an input time series x(t), 
the output time series y{t) will be perfectly correlated to 
x(t) provided that no extra noise enters in the measure- 
ments. In this case the intrinsic coherence function has 
a value of 1. Two completely unrelated time series have 
instead an intrinsic coherence value of 0. The intrinsic 
coherence function is obtained by removing the Poisso- 
nian noise level n from each power spectrum and corre- 
lation function. We then define the intrinsic coherence 
by following Scaringi et al. (2013): 



iLtiv) 



l((CH)| 2 -" 2 
({\X{u)\*)){{\Y{v)\*)y 



(6) 



where the double brackets refer to averages over the N 
ensembles and the number of adjacent Fourier frequency 
bins. 

3. RESULTS 

We detect the 1 Hz modulation in seven distinct Obs- 
Ids and six different outbursts (30-Aug-2009, 21-Mar- 
2010, 12-Jun-2010, 04-Oct-2010, 23-Jan-2011, 06-Nov- 
2011; see Figure 1 and Table 1), with power in a nar- 
row range of spectral frequencies (~ 0.01 — 10Hz). In 
the remaining four outbursts we cannot place any sig- 
nificant constraint on the presence of the 1 Hz modula- 
tion. Indeed one outburst was observed with the Swift 
X-Ray Telescope (XRT) in photon counting mode (2.5 
s time resolution) with a Nyquist frequency of only 0.2 
Hz. The remaining three outbursts (30-Jul-2009, 01-Oct- 
2009, 28-Oct-2009) were observed within the PCA Galac- 
tic Bulge Scan program (Swank & Markwardt 2001) and 
analyzed in Altamirano et al. (2010). In these latter three 
outbursts the power spectra show no significant feature 
due to the poor counting statistics dominated by photon 
counting noise. 

In some observations the 1 Hz modulation is also 
clearly observed in the time domain (see Figure 2 and 
compare to Figure 4 in Patruno et al. 2009a). The 1 Hz 
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modulation appears as a Gaussian in the power spec- 
trum along with one or more other Gaussians whose 
centroid frequency is harmonically related to the "fun- 
damental" Gaussian centroid frequency (sec Figure 3). 
The quality factor of the Gaussians is rather low, with 
Q = vq/FWHM ~ 1 - 2 and the amplitude is 40-55% 
rms (2-60 keV full energy band) in three observations. 
In the remaining observations the feature turns into a 
broadband-limited noise feature with rms amplitude be- 
tween 20 and 40% rms. In all these observations the 
X-ray flux is always above 3 mCrab, whereas in all the 
others with a lower flux, the source NGC 6440 X-2 is 
not detected with a measured flux consistent with back- 
ground emission. 

In addition to the 1 Hz modulation, in two observa- 
tions we also detect a broad QPO (Q ~ 1) with centroid 
frequency of 7.8 ± 1.0 Hz (Obs-Id: 94044-04-02-00; re- 
ported also in Altamirano et al. 2010) and 16.8 ± 0.9 Hz 
(Obs-Id: 94315-01-14-00). This QPO is well fitted in the 
power spectrum with a Lorentzian function. 

3.0.1. Energy Dependence, Coherence and Time Lags 

We detect a clear increase of the 1 Hz modulation rms 
amplitude at high energies, with a rms amplitude that 
grows linearly from 40% up to a peak of ~ 60% rms 
in the 16-28 keV energy band (Figure 4). The intrinsic 
coherence between two broad energy bands (2.47—7.35 
keV and 7.76—28.67 keV) is « 1 at frequencies between 
0.5 and 1.5 Hz, and drops substantially outside this fre- 
quency range. Above 3 Hz and below 0.1 Hz the statistics 
become poor and the data points in the intrinsic coher- 
ence function have large error bars (Figure 5). We use 
the cross-spectra to inspect for the presence of time lags 
between the soft and hard energy bands. The time lags 
are defined as: 



t{v) = 



${v) arg[((C(u)}}} 



2-kv 



2-kv 



(7) 



where ( E>(^) is the Fourier phase lag between two identical 
Fourier frequencies [y) measured between the soft and 
hard energy bands. No time lags are detected between 
the two time series, with typical upper limits on the time 
lags between a specific Fourier frequency of less than 5 
ms. 

3.1. Pulsations 

We detect significant pulsations at a frequency of 
205.892 Hz in all seven Obs-Ids. The frequency shows 
a Dopplcr modulation consistent with the orbital solu- 
tion reported in Altamirano et al. (2010), although we 
do not attempt to refine the solution with a coherent 
timing analysis. The pulsations are detected in coinci- 
dence with the presence of the 1 Hz modulation. 

4. DISCUSSION 

The spectral and timing properties of NGC 6440 X- 
2 provide a wealth of observations to compare with the 
predictions of the trapped disk and disk instability model 
presented in D'Angelo & Spruit (2010, 2011, 2012). Be- 
low we argue that the 1 Hz modulation spectral energy 
dependence and coherence between energy bands suggest 
the variability is driven by a changing accretion rate onto 
the star, while the 1 Hz modulation timescalc, amplitude 
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Fig. 3. — Power spectrum of the observation 94315-01-14-00 in 
the broad energy band (top panel), soft band (middle panel) and 
hard band (bottom panel). The 1 Hz modulation appears in the 
RS 0.1 — 10 Hz band with a peak at about 1 Hz and can be modeled 
with a Gaussian plus several harmonics. A higher frequency QPO 
appears at Ri 16 Hz and is seen with more evidence in the soft 
band, whereas it turns into a broad feature at higher energies. 

and mean accretion rate are in line with the predictions 
from D'Angelo & Spruit (2010, 2012). We also discuss 
the detailed timing properties of the source: the unusual 
power spectral shape of the modulation, and the short 
recurrent outbursts. 

4.1. Evidence for accretion rate modulation 

Our analysis of the outbursts of NGC 6440 X-2 shows 
that the 1 Hz modulation is present with high amplitude 
in at least six out of ten outbursts and that it is a recur- 
rent phenomenon. The lack of the 1 Hz modulation in 
the remaining four outbursts is unconstrained due to the 
poor photon flux or time resolution limits. The timing 
and energy properties of the outburst strongly suggest 
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Fig. 4. — Energy dependence of the 1 Hz modulation fractional 
rms amplitude (Obs-Id: 94315-01-14-00). The fractional amplitude 
increases with energy from 40% rms up to 60% in the higher energy 
band. The errorbars on last point (16-28 keV) are very large but 
are consistent with a linear increase in the rms amplitude of the 
1 Hz modulation. Energies above w 30 keV are not included as 
they are dominated by background noise. The energy dependence 
closely resembles what observed in SAX J1808. 4-3658 (see Patruno 
et al. 2009a). 
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Fig. 5. — Intrinsic coherence function between energy bands 
2.47-7.35 keV and 7.76-28.67 keV for the observation 94315-01- 
14-00 (the observation with the highest signal-to- noise ratio). The 
1 Hz modulation appears in the ps 0.1 — 10 Hz band with a peak 
at about 1 Hz. The intrinsic coherence function shows how the 1 
Hz modulation appearing in the two X-ray lightcurves in the two 
energy bands selected is causally connected (intrinsic coherence 
«1). 

that the modulation is produced by variations in the ac- 
cretion rate onto the star. 

Coherent pulsations at 206 Hz are clearly detected dur- 
ing all outbursts, demonstrating that the channeled ac- 
cretion flow persists when the 1 Hz modulation is seen. 
A complete study of the periodic variability of the pul- 
sations will be presented elsewhere. 

The energy spectrum (Heinke et al. 2010) and energy- 
timing information (this work) also suggests a large- 
amplitude variation in accretion rate onto the star at 
one second intervals. The spectral energy distribution 
across all outbursts can be fit with a power law with an 
index between fa 1.7-3.4 (Heinke et al. 2010), which 



is most likely produced by inverse Compton scattering of 
thermal seed photons, both from hot accretion-shocked 
gas on the stellar surface and hot optically thin gas in 
the inner accretion flow. 

The soft and high energy bands are strongly coherent 
between 0.5-1.5 Hz, consistent with an overall change in 
luminosity produced by a changing accretion rate onto 
the star, where the variability timescale is long enough 
for the optically thin corona above the neutron star and 
thermal emission from the neutron star surface to reach 
quasi-equilibrium (Haardt & Maraschi 1991). The drop 
in coherence above this frequency may then be a reflec- 
tion of the intrinsic timescale for density or temperature 
fluctuations in the hot corona. 

We also measure a steep increase (by about 20%) of 
the 1 Hz modulation rms amplitude with energy, so that 
the flux at higher energies varies more than the flux at 
low energies. Assuming that the overall emission is dom- 
inated by the accretion shock right above the neutron 
star's surface (see e.g., Gierlihski k Poutanen 2005), the 
increased rms amplitude at higher energies naturally re- 
sults from a changing accretion rate. 

Changes in the density and temperature in an opti- 
cally thin corona above the neutron star will lead to a 
variable slope of the power-law index in the spectrum. 
(The slope is set by the Compton y parameter, a func- 
tion of the gas's temperature and optical depth, where 
higher temperatures and higher densities both produce 
harder spectra) . A fluctuation in accretion rate will thus 
change not only the overall flux but also the hardness 
of the inverse Compton spectrum, so that flux at higher 
energies is more variable than at lower energies, as seen 
in both NGC 6440 X-2 and SAX J1808.4-3658. 

No thermal component was detected in the 49. Iks 
Chandra observation during the July 28, 2009 outburst 
(Heinke et al. 2010) or any other (less constraining) ob- 
servation. This is in contrast to SAX J1808. 4-3658, in 
which a thermal component is frequently seen in out- 
bursts, and is attributed to a hot spot on the surface of 
the neutron star (Poutanen k Gierlihski 2003; Gierlihski 
k Poutanen 2005; Patruno et al. 2009b; Papitto et al. 
2009; Kajava et al. 2011). In SAX J1808.4-3658 the 
thermal component is observed to lag the high energy 
one by ~ 0.125 ms (Cui et al. 1998; Gierlihski et al. 
2002), which Gierlihski k Poutanen (2005) attribute to 
strong Doppler beaming of the blackbody emission com- 
pared to the more diffuse inverse Compton emission. The 
lack of thermal emission and low pulsation fraction may 
indicate that the main accretion hot spot in NGC 6440 
X-2 is beamed out of our line of sight, while the upper 
limits on phase lags between hard and soft energy bands 
(of the order of 5 ms) are not particularly constraining. 



4.2. Oscillations ir 



unstable dead disk 



The instability studied in Spruit k Taam (1993) and 
D'Angelo k Spruit (2010, 2012) appears when the in- 
ner edge of the disk moves into a marginally dead state, 
with r m ~ r co . This condition requires a low mean ac- 
cretion rate so that the inner edge of the disk moves 
temporarily outside the corotation radius. NGC 6440 
X-2, which has one of the lowest outburst luminosities 
of all known AMXPs, supports this picture. Its typi- 
cal X-ray outburst luminosity is « 10 36 ergs -1 (Heinke 
et al. 2010), corresponding to a mass accretion rate of 
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w 2 x lO~ lo M0yr -1 for a canonical neutron star mass 
of 1.4 Mq and radius of 10 km. 

The disk truncation radius r m is set by the magnetic 
field, and is subject to considerable uncertainty due to 
the uncertainty in exactly how the disk and field interact. 
We follow Spruit & Taam (1993) and D'Angelo & Spruit 
(2010) and define the inner disk radius as where the field 
is strong enough to force the gas to co-rotate with the 
star: 



r m = 47.7k 




B„ 



2/5 



M* 



( — ) 

V 10km J \ 2 x lO-^il/oyr- 1 



4xl0 8 G/ V!- 4M 

\ -1/5 

M 



-1/10 

(8) 



-V5 „ x 3/10 



4.9ms 



where the subscript "*" identifies quantities referred to 
the neutron star and B^jB* is the fractional size of the 
toroidal field induced by the relative rotation between 
the magnetic field and the disk. For B^/B* > 1, the 
field lines will rapidly inflate and sever the connection 
between the disk and the star (e.g., Aly 1985; Lovelace 
et al. 1995), but the exact average value is uncertain. 
Assuming B<f,/B* ~ 0.1 — 1 implies an uncertainty of 
up to ~ 60% on the value of -B* that will put r m ~ r c . 
The inferred magnetic field for NGC 6440 X-2 (assuming 
B^/B* = 0.1 and r m = r co for M ~ lO" 10 M yr- 1 ) will 
be 4 x 10 8 G, well in line with most AMXPs (see e.g., 
Patruno & Watts 2012 for a review). 

The timescale of the 1 Hz modulation also suggests an 
instability in the inner accretion disk. This timescale is 
much longer than the spin frequency of the star (206 Hz) 
or the dynamical timescale of the inner disk (similar 
when r m ~ r co ). It is also likely much longer than the 
timescale on which the field lines are twisted open and re- 
connect to the disk, although this is currently somewhat 
uncertain (e.g., Romanova et al. 2009). 

The disk instability studied by Spruit & Taam (1993) 
and D'Angelo & Spruit (2010, 2012) causes density fluc- 
tuations in the inner parts of the disk, which generally 
evolve on viscous timescales. The instability frequency 
roughly scales with the viscous timescale around the in- 
ner edge of the disk. This timescale is given by: 



T v i S c ~ 80 ( — 
a 



-4/5 



M 



2 x 10" 10 Af Q yr- 



-3/10 



M* 



1.4 M P 



© 



1/4 



50 km 



5/4 



(9) 



where a ~ 0.1 is the canonical a-disk parameter for 
geometrically thin and optically thick disks (Siuniacv 
& Shakura 1977) and r ~ r m . For NGC 6440 X-2, 
Tvisc ~ 80 s, and the modulation at 1 Hz is well in the in- 
stability range seen by D'Angelo & Spruit (2010, 2012). 
In section 4.4, we discuss this further and place tighter 
physical constraints on the disk-magnetosphere interac- 
tion based on the instability properties. 

4.3. A dead disk in quiescence? 

As well as being weak with short durations, the out- 
bursts seen in NGC 6440 X-2 have the shortest recur- 
rence time of any known AMXP, with ten outbursts 



seen at roughly monthly intervals between 2009 and 
2011. Similar behaviour was seen in another AMXP, 
IGR J00291+5934, which in 2008 underwent two small 
outbursts (lasting around 10 days each) separated by 
a month in quiescence (Patruno 2010; Hartman ct al. 
2011). This AMXP also showed a QPO at 0.5 Hz (Hart- 
man et al. 2011), although this QPO is much more typical 
(with substantially lower amplitude of (13 ± 2)%, and 
a standard Lorentzian shape) than those seen in NGC 
6440 X-2 and SAX J1808.4-3658. There also may be 
a connection between the appearance of the instability 
in SAX J1808.4-3658 and the 'reflaring' portion of the 
light-curve. The 'reflaring' refers to large luminosity vari- 
ations (up to 3 orders of magnitude) on timescales of 3-5 
days in the decay tail of SAX J1808.4-3658 (see e.g., 
Figure 1 in Patruno et al. 2009a), suggesting large-scale 
changes in the accretion rate as the source approaches 
quiescence. 

The short variation timescale of the reflares and out- 
bursts is incompatible with the standard ionization in- 
stability model believed to power normal outbursts (e.g., 
Lasota 2001) and suggests that a substantial amount 
of gas remains in the inner regions of the accretion 
disk (close to r co ) when the accretion rate decreases, 
exactly as is expected for a dead disk. As was first 
noted by D'Angelo & Spruit (2012), the unusual fre- 
quent outbursts and strong 1 Hz modulation are both 
well-explained by a disk that never moves into a strong 
propeller phase at very low M, but instead stores mass 
as a 'dead disk', which becomes unstable at higher M. 

The density structure of a dead disk, with its low tem- 
perature and a substantial amount of stored mass in 
the disk inner regions, is considerably different from the 
standard accretion disk and could strongly alter the ion- 
ization instability cycle. Although detailed calculations 
of how the ionization instability would be triggered in 
a dead disk are beyond the scope of the current work, 
the additional stored mass during quiescence could lower 
the density threshold for which newly accumulated mass 
(from the star or the outermost regions of the accretion 
disk) can trigger a new outburst. Such outbursts would 
be expected to have both shorter recurrence times and 
weaker outbursts, since less additional mass is accumu- 
lated between outbursts. 

A further prediction of the dead disk picture is that 
spin-down should be enhanced during the quiescent 
phase, since the disk inner edge remains close to r co and 
continues to spin down the star. Note however, that the 
disk-field coupling parameters we infer for this source (in 
section 4.4) imply a narrow connected region between 
the star and the disk, which will also limit the spin-down 
from this source. 

4.4. Constraining properties of the disk-field interaction 

In the model presented in D'Angelo & Spruit (2010, 
2012), the instability develops in two different regimes, 
depending mainly on the mean mass accretion rate 
through the disk. The more-studied regime (called RI 
in D'Angelo & Spruit 2012) occurs at very low accre- 
tion rates, with M w 10~ 4 - 0.1M CO (where M co is the 



accretion rate that sets r c 



It is characterized 



by long periods of quiescence while mass accumulates in 
the inner regions of the disk followed by brief accretion 
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Fig. 6. — Illustration of the interaction region Ar where the 
neutron star and disk magnetic fields are coupled. At r < r m the 
neutron star magnetic field lines arc closed, whereas at r m < r < 
r m +Ar the field lines open and interact with the accretion disk. 
Beyond the interaction region the disk flow is unperturbed. 




Fig. 7. — Transition region Ar2 around r co . When r m is smaller 
than r co — — j 3 - the gas accretes along magnetic field lines, and 
spins up the star. When r m lies in the transition region r co ± J 2 
the disk is in a marginal state, where accretion onto the star con- 
tinues but the inner disk simultaneously begins extracting angular 
momentum from the star, so that the spin up turns to spin down. 

When r m lies above r co H ^2. the accretion flow stops but the 

inner disk radius remains close to the corotation radius. When 
the 1 Hz modulation is seen, the inner disk edge stays around 

r ~ r + Ar 2 
' m — 'co ^ o 



outbursts. Indeed, the authors found that the period of 
the instability could span several orders of magnitude, 
from around 0.1 — 10 T v i sc (where T v i sc is the viscous 
timescale at r m , Equation 9), depending chiefly on the 
time needed to fill up the reservoir and tip the source 
into an accretion burst. 

The instability outburst shape and duration also de- 
pends quite sensitively on the details of the disk-field 
interaction, which D'Angelo & Spruit (2010, 2012) 
parametrized as two length-scales: the width of the cou- 
pled region between the field and the disk, Ar, and Ar2, 
the transition length around r co between the accret- 
ing and non- accreting solutions (see Figure 6 and 7). 
Both of these length scales are assumed to be small: 
Ar/r co , Ar2/r co <C 1. In the RI instability region, 
D'Angelo & Spruit (2010, 2012) found that the insta- 
bility occurs for a broad range of Ar/r m (see Figure 3 
and 4 in D'Angelo & Spruit 2012) but only when Ar2/r m 
is very small (of order of 0.01). 

However, in D'Angelo & Spruit (2012) the authors 
identified a second instability region, which they referred 



to as RII. Here the mean accretion rate is 0.1 — 10 x M co - 
much higher than RI - and accretion continues through- 
out the cycle, so that the amplitude of the instability 
is much smaller. The instability timescale is also much 
shorter than for RI, with periods around 0.01 — 0.1 T v i sc , 
and has a more sinusoidal, regular shape and a smaller 
amplitude, although the outburst length is still typi- 
cally shorter than the build-up time (see e.g., Fig. 5 
of D'Angelo & Spruit 2012). The instability also ap- 
pears for different disk-field interaction parameters: here 
Ar/r m is very small (« 0.05), while Ar^/rm has a much 
broader range of instability. 

In NGC 6440 X-2 the corotation radius is located at 
approximately 50 km from the neutron star center, but 
the surface magnetic field is still unknown, so that it 
is not possible to estimate the location of the magne- 
tosphcric radius from the luminosity alone. However, 
if we assume a typical AMXP magnetic field strength 
of ps (1 — 5) x 10 8 G, the observed mass accretion rate 
is of the same order of magnitude or lower than M co , 
suggesting that we are observing the RII instability. In 
SAX J1808.4-3658 (with B w 10 8 G, see Hartman et al. 
2008; Patruno ct al. 2012), the 1 Hz modulation was 
seen at fluxes corresponding to M/M co w 1, suggesting 
that source is also unstable in the RII region. In this 
regime the instability only appears over a relatively nar- 
row range of M, offering a plausible explanation for why 
the instability is only observed in the decay tail of the 
outburst. 

One obvious point of interest is why the 1 Hz modu- 
lation central frequency is nearly the same in both SAX 
J1808.4-3658 and NGC 6440 X-2, given that the system 
parameters are quite different. Despite the comparable 
X-ray flux detected in both AMXPs when the 1 Hz modu- 
lation is observed, the corresponding mass accretion rate 
of SAX J1808. 4-3658 is about an order of magnitude 
lower than in NGC 6440 X-2 because the first source lies 
at about 3.5 kpc (Galloway & Cumming 2006) whereas 
the latter is at 8.5 kpc (Ortolani et al. 1994). If we 
compare the fundamental centroid frequency of the 1 Hz 
modulation observed in all observations of NGC 6440 X- 
2 with those of SAX J1808. 4-3658 then it is, on average, 
systematically lower in NGC 6440 X-2 by only a factor 
of pa 2. Without knowing the magnetic field in cither 
source it is not easy to constrain exactly where the disk 
is truncated (or, to invert the problem, the actual ratio 
M/M co ), but if we assume that the relevant timescale is 
the viscous timescale at r co then this will be ~ 80 s for 
NGC 6440 X-2 vs. ~ 40 s for SAX J1808.4-3658, so the 
model predicts roughly the same v for both sources, as 
observed. An independent constraint on B* would allow 
a much more secure determination of M/M co and also 
Ar2 and Ar. 

The RII instability appears only for very small values 
of Ar, of order 0.01-0.05 r m . This suggests a very small 
coupling region between the field and the disk - of order 
1 km - in both SAX J1808.4-3658 and NGC 6440 X- 
2. Given the uncertainties in determining the system's 
physical parameters, Ar2 cannot be as well-constrained, 
although at these accretion rates Ar2 cc j/ _1 (where v 
is the instability frequency; see Fig. 7 of D'Angelo & 
Spruit 2012), so a smaller value for Ar2 (of about 1-10 
km) would also seem to be favored. 
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D'Angclo & Spruit (2012) found that when the insta- 
bility is triggered in the RII regime, the accretion torques 
are reduced compared with the standard disk solution. 
An implication of this is that if an AMXP is observed 
to show the 1 Hz modulation in a certain mass accretion 
rate range, then the observed accretion torques must be 
smaller than the prediction made with the approxima- 
tion: 

N = My/GM*r m (10) 

It is possible to speculate therefore that the lack of strong 
positive accretion torque in SAX J1808. 4-3658 (Hart- 
man et al. 2008, 2009) might be partially caused by the 
source being almost always close to this marginally ac- 
creting state. The fact that SAX J1808.4-3658 is found 
with r m close to r co was also suggested by Haskell & 
Patruno (2011) based on coherent timing and quiescent 
thermal properties of the accreting neutron star. 

4.5. Detailed timing properties 

The fact that the 1 Hz modulation can be fitted with a 
Gaussian and not with a Lorentzian suggests some intrin- 
sic difference between the generation of the 1 Hz modu- 
lation and the higher frequency QPOs observed in NGC 
6440 X-2 and SAX J1808.4-3658. It also suggests some 
remarkable difference between most of the other QPOs 
broadly observed in other low mass X-ray binaries and 
the 1 Hz mechanism. 

If a Gaussian stochastic process controls the fluctua- 
tions in mass accretion rate around the magnctosphcric 
radius, then the expected M profile instability will be 
a Gaussian. This in turn might explain why the 1 Hz 
modulation has a Gaussian shape and not a Lorentzian. 
It is conceivable indeed to think that, once the insta- 
bility is triggered and the gas is accumulating around 
I'm ~ r co , the main timescale of the instability will be 
strictly periodic if none of the quantities in the system is 
changing (i.e., M, Ar/r m and Ar2/V m ). However, if we 
let M fluctuate (in a Gaussian stochastic way) around its 
mean value, the periodicity of the instability will slightly 
change, thus broadening the periodicity across adjacent 
frequencies. It is then possible to understand why the 
shape of the 1 Hz modulation is Gaussian and why it 
becomes a broad feature in some observations whereas 
in others it is a more periodic and narrow feature: the 
stronger the fluctuations in M around the mean value, 
the broader the 1 Hz modulation will appear in the power 
spectra. 

4.6. Alternative Models 

Patruno et al. (2009a) reported an extensive discus- 
sion of alternative models to the trapped-disk instabil- 
ity presented in Spruit & Taam (1993) (and later devel- 
oped in D'Angelo k Spruit 2010, 2012). It was shown 
that none of the alternative models can explain at the 
same time the observed features of the 1 Hz modula- 
tion in SAX J1808. 4-3658. Identical consideration apply 
also to the 1 Hz modulation in NGC 6440 X-2, as all 



the observational features closely resemble those in SAX 
J1808. 4-3658 and identical considerations apply. More 
recently, Romanova ct al. (2013) has carried 3D MHD 
simulations with a dipolc magnetic field misaligned from 
the neutron star's rotational axis, in which a warped 
disk appears. Perturbations of the disk in the vertical 
direction produce low-frequency bending waves, whose 
frequency is approximately equal to the Keplerian fre- 
quency (Cl = ^GAI/r 3 1 of the gas in the disk. Since 

the inner disk is truncated at about ~ 50 km, the Ke- 
plerian frequency there is of the order of 200 Hz, which 
is 2-3 orders of magnitude off the observed instability 
frequency. The warped disk predicts oscillations with 
frequency down to 0.1-1 Hz, but the region of the disk 
that oscillates has to be far away from the neutron star, 
at several thousands kilometers, where no X-ray radia- 
tion is produced and the inner accretion flow (at about 
50 km) would still be modulated at the much higher Ke- 
plerian frequency there. 

5. CONCLUSIONS 

We have studied the 1 Hz modulation in the AMXP 
NGC 6440 X-2 and suggested that this feature has the 
same physical origin as seen in the other AMXP SAX 
J1808. 4-3658. We propose the identification of this fea- 
ture with a magnctosphcric instability triggered in the 
weak propeller stage. Our study suggests a very small 
transition region Ar/r m of the order of 0.01-0.05. This 
region is where the magnetic field lines open and inter- 
act with the accretion disk and this has to be small in 
both AMXPs considered in this work. We are not able to 
put such strong constraints on transition region Ar2/r m , 
where the nature of the accretion torques changes from 
spin up to spin down, but suggest it lies in the range 0.01- 
1. Given that r m ~ r co when the instability is observed, 
the radial extension of the two transition regions for the 
two AMXPs is: Ar ps 1 km and Ar 2 ~ 1 - 10 km. This is 
the first time that the width of the magnctosphcre-disk 
interaction region can be constrained with direct obser- 
vations. 

We have also suggested an explanation for the Gaus- 
sian shape of the 1 Hz modulation as seen in the power 
spectra as due to fluctuations in the mean mass accretion 
rate during the instability. Strong M fluctuations can 
smear the 1 Hz modulation from a rather coherent QPO 
to a broadband-limited noise. If confirmed, this might 
be the first clear physical identification of the mecha- 
nism behind the generation of a low frequency QPO and 
broad band noise. 
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